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I  Introduction  Distribution  Unlimited 

This  research  project  focused  on  issues  met  in  using  large-eddy  simulation  (LES)  to  predict 
turbulent  nonpremixed  combustion.  LES  has  been  recognized  as  a  very  promising  approach  to 
modeling  such  flows[l,  2,  3],  In  this  approach  the  subgrid-scale  (SGS)  scalar  mixing  and  the 
resulting  instantaneous  distribution  of  scalar  values  in  each  grid  volume  (i.e.,  the  filtered  density 
function)  must  be  faithfully  represented  in  order  to  accurately  predict  the  chemical  reaction  rate. 

An  important  modeling  method  uses  the  transport  equation  of  the  filtered  joint  density  function 
of  velocity  and  scalars  in  which  both  the  reaction  source  term  and  the  turbulent  transport  terms 
are  in  closed  form[3].  This  research  studied  issues  in  using  this  approach  by  investigating  the  SGS 
mixing  of  conserved  scalars  which  often  play  a  crucial  role  in  LES  of  nonpremixed  combustion. 

Significant  progress  has  been  made  in  understanding  the  SGS  velocity  and  scalar  filtered  density 
function  and  their  dynamics  and  issues  of  LES  modeling.  Specifically,  the  following  work  was 
completed: 


•  Proved  Galilean  invariance  of  the  velocity  probability  density  function  (PDF),  velocity  fil¬ 
tered  density  function,  their  transport  equations,  and  the  statistics  conditional  on  one-point 
velocity.  The  proof  provides  a  formal  justification  of  using  the  PDF-based  approach  to  model 
turbulent  flows  in  a  manner  that  is  invariant  in  any  inertial  reference  frame. 
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•  Investigated  the  filtered  joint  density  function  of  velocity  components  (VFJDF)  and  its  trans¬ 
port  equation  using  experimental  data.  The  VFJDF  and  the  structure  of  the  SGS  velocity 
are  found  to  strongly  depend  on  to  degree  of  nonequilibrium  spectral  transfer. 

•  Investigated  the  dependence  of  the  structure  of  SGS  velocity  on  the  nonequilibrium  spectral 
energy  transfer.  The  results  are  important  for  establishing  a  framework  for  modeling  the 
degree  of  nonequilibrium  spectral  transfer  and  the  VFJDF. 

•  Investigated  the  filtered  joint  density  function  of  velocity  components  and  a  conserved  scalar 
(VSFJDF)  and  its  transport  equation  using  experimental  data.  The  VSFJDF  and  the  cou¬ 
pling  between  the  SGS  velocity  and  scalar  are  found  to  depend  strongly  on  the  degree  of 
nonequilibrium  spectral  transfer. 

•  Investigated  the  dependence  of  the  structure  of  SGS  scalar  fields  on  nonequilibrium  spectral 
energy  and  scalar  transfer.  The  dependence  is  important  for  establishing  a  framework  for 
modeling  the  degree  of  nonequilibrium  spectral  transfer  and  the  VSFJDF. 

•  Investigated  the  issues  of  modeling  the  filtered  energy  dissipation  and  the  filtered  scalar 
dissipation. 

•  Conducted  a  preliminary  study  of  SGS  mixing  of  the  mixture  fraction  and  the  SGS  flame 
structure.  The  study  was  based  on  the  knowledge  gained  in  the  above-mentioned  work,  and 
the  preliminary  results  have  strong  implications  for  studying  SGS  flame  structure  and  local 
extinction  and  reignition. 

The  primarily  findings  of  this  research  are  discussed  in  the  following. 
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(joint)  probability  density  function 
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SGS  scalar  variance  • 
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Ill  Findings 


The  velocity-scalar  FJDF  method  was  recently  established  by  Peyman  Givi  and  Stephen  Pope 
[3].  It  solves  the  velocity-scalar  FJDF  transport  equation  in  which  the  SGS  scalar  transport  term 
is  in  closed  form.  More  importantly,  this  approach  allows  mixing  models  to  depend  on  velocity, 
thus  having  the  potential  to  provide  a  more  realistic  description  of  the  SGS  mixing,  which  is  local 
in  physical,  velocity,  and  scalar  spaces.  Consistency  with  these  conditions  of  localness  is  key  to  the 
accuracy  of  PDF-based  approaches.  We  studied  the  FJDF  and  its  transport  equation  for  improving 
the  velocity-scalar  FJDF  method. 


The  FJDF  of  velocity,  u,  and  scalar,  <f>  is  defined  as 
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fu<t>{y,  <f>] x,t)  =  J  n<Sh(x',t)  -  Vi]S[(j>(-x!  ,t)  -  4>\G(x'  -  x)dx', 
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where  t,  x,  v,  <j>,  5 ,  and  G  are  time,  the  physical-space  coordinates,  the  sample-space  variables  for 
u  and  <f>,  the  Dirac  delta  function,  and  the  filter  function,  respectively.  The  integration  is  over  all 
physical  space.  The  FJDF  represents  the  weighted  joint  distribution  of  the  velocity  components  and 
the  scalar  in  a  grid  cell.  The  transport  equation  of  the  FJDF  is  obtained  using  the  Navier-Stokes 
equations  and  the  scalar  transport  equation: 
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where  (  •  |u  =  v,  <j)  =  4>)i,  denotes  a  conditionally  filtered  variable  conditional  on  the  velocity 
vector  and  the  scalar.  The  left  hand  side  represents  the  time  rate  of  change  of  the  FJDF  and 
transport  of  FJDF  in  physical  space.  The  terms  on  the  right  hand  side  are  transport  of  the  FJDF  in 
velocity  space  by  the  resolvable-scale  pressure  gradient,  SGS  pressure  gradient,  viscous  acceleration, 
and  molecular  diffusion.  An  alternative  term  to  the  conditionally  filtered  scalar  diffusion  is  the 
conditionally  filtered  scalar  dissipation  (x|u  =  v,  (j) 
dissipation  rate. 


,  r^dd>  dd>  . 
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Unlike  a  PDF,  which  is  a  statistic,  the  FJDF  is  a  random  process,  and  therefore  requires  statisti¬ 
cal  descriptions.  In  this  research  we  employed  conditional  sampling  and  averaging  techniques.  Our 
previous  studies  of  conserved  scalar  FDF  and  the  FDF  transport  equation  using  these  techniques 
[4,  5,  6]  have  shown  that  the  SGS  scalar  has  close-to-Gaussian  and  bimodal  FDFs  for  small  and 
large  instantaneous  SGS  variance,  respectively.  In  addition  the  SGS  scalar  is  in  spectral  equilib¬ 
rium  (production  of  SGS  variance  smaller  than  dissipation)  and  non-equilibrium  (production  larger 
than  dissipation),  respectively.  The  bimodal  SGS  scalar  contains  diffusion  layers,  thus  is  similar 
to  a  nonpremixed  scalar  field.  The  results  suggest  that  mixing  models  that  can  better  predict  the 
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non-equilibrium  regime  are  needed.  In  this  work  we  investigated  the  statistical  dependence  between 
the  SGS  velocity  and  scalar  fields,  and  the  effects  of  the  SGS  velocity  on  SGS  scalar  mixing.  The 
FJDF  and  the  SGS  terms  in  the  FJDF  transport  equation  are  analyzed  using  their  conditional 
means  with  the  SGS  scalar  variance  (</>/,2)l  and  the  SGS  kinetic  energy,  kL  =  ((u"2)L  +  (u'^)L)/2 
as  conditioning  variables. 

We  conducted  experiments  in  an  axisymmetric  turbulent  jet  and  used  the  experimental  data  to 
analyze  the  unclosed  terms  in  the  velocity-scalar  FJDF  and  its  transport  equation.  Measurements 
were  made  at  80  jet  nozzle  diameters,  D,  (=15  mm)  downstream.  The  jet  Reynolds  number  UD/v 
was  40000.  A  sensor  array  consisting  of  three  hot-wire  probes  was  used  to  obtain  data  for  performing 
filtering  operations  in  both  the  streamwise  (a:)  and  transverse  (y)  directions.  The  filter  size  was 
varied  from  63  to  250  Kolmogorov  scales.  This  technique  enabled  us  to  acquire  a  large  amount  of 
data  (2  x  108  samples)  necessary  to  achieve  sufficient  statistical  convergence. 

In  the  following  we  briefly  outline  the  results  from  our  study  and  their  implications  for  LES. 
A.  Galilean  invariance  of  velocity  PDF  and  FJDF 

In  turbulence  and  turbulent  combustion  modeling  it  is  generally  necessary  for  the  model  equa¬ 
tions  to  have  the  same  Galilean  transformation  properties  as  the  exact  transport  equations.  The 
Navier-Stoke  equations  and  the  moment  equations  for  turbulent  fluctuations  are  known  to  be 
Galilean  invariant.  Thus  the  modeled  equations  must  also  be  Galilean  invariant. 

PDF  methods  have  become  an  important  modeling  approach  in  turbulent  flows,  especially 
for  reacting  flows[7,  1].  In  this  approach  the  PDF  transport  equations  are  modeled.  Although  the 
modeled  velocity  PDF  equations  are  implicitly  assumed  to  be  Galilean  invariant  in  the  literature [8], 
the  Galilean  transformation  properties  of  both  the  exact  and  modeled  PDF  equations  have  not  been 
given  previously  (even  the  transformation  property  of  the  PDF  itself  has  not  been  discussed). 

Another  motivation  for  this  analysis  was  the  Galilean  transformation  properties  of  statistics 
conditional  on  turbulent  velocity,  which  are  important  in  studying  turbulent  flows.  In  fact,  the 
key  terms  in  the  one-point  PDF  equation  are  means  conditional  on  velocity  at  a  point.  Statistics 
conditional  on  the  resolvable-scale  velocity  are  also  frequently  used  in  analyses  of  the  subgrid-scale 
stress  and  models  in  large  eddy  simulation.  However,  the  Galilean  invariance  issue  associated  with 
such  conditional  statistics  has  not  been  previously  addressed  and  there  seems  to  be  some  confusion 
in  the  literature  as  to  whether  or  not  they  are  invariant  under  Galilean  transformations.  In  this 
work  we  also  considered  this  issue. 

We  first  examined  the  transformation  properties  of  the  velocity  joint  probability  density  function 
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(JPDF) 
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/«(v;  x,  t)  =  (J]  5 [tii  -  wj)  (3) 

1=1 

where  an  (•)  denotes  an  ensemble  mean.  The  Galilean  transformations  from  one  inertial  reference 
frame  (x,  t)  to  another  (x*,  t*)  are 

x*  =  x  +  Ut,  t*  =  t  (4) 

where  U  is  the  (constant)  velocity  of  the  reference  frame  (x,t)  relative  to  (x*,t*).  Thus  a  velocity 
field  transforms  as 

u*(x*,  t *)  =  u*(x  +  Ut,  t)  =  u(x,  t)  +  U  (5) 

Thus  the  velocity  is  not  invariant  under  Galilean  transformations. 

To  transform  the  JPDF,  in  addition  to  the  transformations  in  Eqs.  4  and  5,  a  transformation 
of  the  sample  space  is  needed  because  the  JPDF  is  also  a  function  of  v.  Since  an  event  u  =  v  is 
equivalent  to  u*  =  v  +  U  in  the  new  frame,  the  transformation  is 

v*  =  v  +  U.  (6) 

Thus  the  sample  space  is  not  invariant.  The  transformation  of  the  JPDF  then  is 

3  3  3 

/.(v; X. ()  =  <n  «fm  -  «i]>  =  <n  <[«?  -vt-  w  -  !/()]>  =  <n  «K  -  «*]}  =  /„•  (v>; **,**)■  P) 

i=l  i— 1  i=l 

3 

Therefore,  the  JPDF  (as  well  as  the  fine-grained  JPDF  JJ  —  Uj])  is  invariant  under  Galilean 

i=i 

transformations  even  though  the  velocity  is  not.  The  invariance  of  the  JPDF  is  due  to  the  invariance 
of  Ui  —  Vi  in  which  the  non-invariance  of  iq  and  Vi  is  canceled.  Since  the  transformations  of  both 
the  velocity  and  the  sample  space  are  involved,  the  transformation  of  the  JPDF  is  conceptually 
somewhat  different  from  those  we  usually  consider,  such  as  the  Reynolds  stress  and  subgrid-scale 
stress.  An  intuitive  way  to  understand  the  invariance  of  the  JPDF  is  that  one  obtains  the  same 
JPDF  values  using  the  variables  before  and  after  the  transformations.  Equation  7  shows  that  the 
functional  form  of  fu»  is  the  same  as  fu  shifted  by  the  amounts  U  and  Ut  in  the  velocity  sample 
space  and  in  physical  space,  respectively. 

Because  the  sample-space  variable  is  an  independent  variable,  the  transformation  is  in  a  way 
similar  to  that  for  a  scalar  variable  <f>.  We  note  that  the  transformation  of  the  sample  space  is  a 
key  step  in  establishing  the  Galilean  transformation  property  of  the  JPDF.  If  this  step  is  omitted, 
the  transformation  group  is  incomplete,  and  a  transformation  v*  =  v  is  implied.  The  incomplete 
group  results  in  the  transformation  /u(v;x, t)  —  fu> (v*  +  U;x*,f*).  Thus  the  JPDF  appears  as 
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non-invariant.  The  transformations  in  Eqs.  6  and  7  preserve  the  transformation  properties  of  the 
moments  of  the  JPDF  (for  example,  the  Reynolds  stress  (u'u'.}  =  (u*'uj')). 


The  JPDF  transport  equation  is  obtained  from  the  Navier-Stokes  equations  and  the  definition 
Eq.  3[7] 

dfu  ,  dfu  d  f/dp,  d  d2Ui 

~9t+ViBxj  =  W"  =  V>'«>  -  (8) 


where  (  •  |u  =  v)  is  a  mean  conditional  on  the  velocity.  Using  the  transformations  in  Eqs.  2,  3,  8 
and  9,  the  transformations  of  the  terms  in  the  JPDF  equation  can  be  obtained  as  follows: 

dfu  dfu .  dfu*  dx*  dfu*  dfu .  dfu  dfu •  (  . 

dt  dt*  dx*,  dt  dt*  *  dx*,  ’  dxj  dx*,  ’  ' 

|-{<^I„  =  v)/4  =  A{(|1|u.  =  v.,/u.}  (10) 

= -k  (11) 

Therefore  the  JPDF  transport  equation  is  invariant  under  Galilean  transformations.  The  invariance 
might  not  seem  intuitive  since  the  conditional  averages  in  Eqs.  10  and  11  shift  in  the  velocity 
sample-space  under  Galilean  transformations.  However,  because  v  is  an  independent  variable, 


dependence  on  it  is  generally  expected  and  does  not  lead  to  non-invariance  as  long  as  the  forms  of 
the  transformed  terms  are  invariant.  Thus  both  the  JPDF  equation  and  its  solutions  (as  well  as  the 
initial/boundary  conditions)  are  Galilean  invariant,  In  contrast,  while  the  Navier-Stokes  equations 
are  Galilean  invariant,  their  solutions  are  not. 


The  conditional  means  in  Eqs.  10  and  11  are  Galilean  invariant  because  the  conditioning  event 
u  =  v  is  equivalent  to  u*  =  v  +  U  =  v*  in  the  new  reference  frame.  This  argument  also  leads  to 
Galilean  invariance  of  the  conditional  statistics  of  any  Galilean  invariant  variable  conditional  on 
the  velocity.  The  invariant  property  provides  an  important  basis  for  using  conditional  statistics 
conditional  on  the  velocity  in  studying  turbulent  flows  and  for  extending  results  to  other  inertial 
reference  frames.  Again,  the  key  step  leading  to  the  invariance  is  the  transformation  of  the  sample 
space.  Like  the  velocity  JPDF,  the  conditional  statistics  appear  as  non-invariant  if  this  transfor¬ 
mation  is  omitted.  This  fact  might  be  a  cause  for  the  confusion  in  the  literature  as  to  whether  the 
conditional  statistics  are  Galilean  invariant  or  not. 


The  above  analyses  were  extended  to  the  filtered  joint  density  function  (FJDF)  used  in  PDF- 
based  large-eddy  simulation  and  its  transport  equation[2).  Speziale[9]  showed  that  the  filtered 
(large-scale)  part  of  a  Galilean  invariant  function  is  also  Galilean  invariant.  Because  the  fine¬ 
grained  JPDF  is  Galilean  invariant,  its  filtered  part,  the  FJDF 


J  ^[ui(x',  t)  -  Vi\G(x.'  —  x)dx'  =  J  <5[u*(x*',  t*)  —  u*]G(x*'  -  x*)dx* 
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is  also  Galilean  invariant.  Thus  all  Galilean-invariant  models  for  the  PDF  methods  can  be  extended 
to  PDF-based  large-eddy  simulation. 

In  summary,  the  JPDF  and  its  transport  equation  were  shown  to  be  Galilean  invariant.  For 
constant  density  flows  they  are  also  invariant  under  extended  Galilean  transformations.  In  addition, 
the  conditional  means  in  the  JPDF  equation  as  well  as  the  conditional  statistics  of  any  Galilean 
invariant  variables  conditional  on  the  velocity  are  invariant.  A  key  step  in  establishing  these 
transformation  properties  is  the  proper  transformation  of  the  sample-space  variable  which  makes 
the  group  of  Galilean  transformations  complete.  The  present  work  provides  support  for  PDF 
methods  and  a  basis  and  justification  for  studying  turbulent  flows  and  turbulent  combustion  using 
statistics  conditional  on  the  velocity. 

B.  Velocity  FJDF,  conditionally  filtered  viscous  acceleration,  and  dissipation 

We  now  discuss  the  results  on  the  velocity  FJDF  and  its  transport  equation.  The  understanding 
gained  provides  a  basis  for  studying  the  velocity-scalar  FJDF. 

Unlike  a  PDF,  a  VFJDF  is  not  a  statistic,  but  a  random  process  and  must  be  characterized 
statistically.  For  a  filter  size  smaller  than  the  integral  length  scale,  the  mean  VFJDF  approximately 
equals  the  velocity  JPDF,  which  is  generally  not  far  from  joint-Gaussian  in  regions  of  fully  devel¬ 
oped  flows  without  large-scale  intermittency[10,  11].  Therefore  the  mean  VFJDF  does  not  provide 
additional  information  on  the  VFJDF.  However,  other  important  characteristics  of  the  VFJDF 
can  be  revealed  by  its  conditional  means.  Recent  results  obtained  by  the  PI  on  conditional  scalar 
FDF  and  scalar-scalar-dissipation  FJDF  have  provided  important  characteristics  of  the  SGS  scalar 
field  [4,  5,  6].  The  SGS  scalar  has  been  found  to  be  in  quasi-equilibrium  (SGS  production  equal  to 
or  smaller  than  dissipation)  and  non-equilibrium  (SGS  production  exceeds  dissipation)  for  small 
and  large  SGS  scalar  variance,  respectively [4,  5].  Here,  for  convenience,  we  refer  to  both  the  cases 
of  SGS  production  equal  to  and  smaller  than  dissipation  as  quasi-equilibrium  because  the  SGS 
scalar  has  very  similar  characteristics  for  the  two  cases.  The  FDF  and  the  terms  in  its  transport 
equation  in  these  two  regimes  have  qualitatively  different  characteristics.  For  equilibrium  SGS 
scalar  fields  the  FDF  is,  on  average,  close  to  Gaussian  (Fig.  1)  and  the  scalar  dissipation  has  only 
moderate  dependence  on  the  SGS  scalar.  These  properties  are  similar  to  those  of  unfiltered  scalar 
fluctuations  in  fully  developed  turbulent  flows;  however,  for  non-equilibrium  SGS  scalar  fields  (at 
large  SGS  scalar  variance),  the  FDF  is  bimodal,  and  the  scalar  dissipation  depends  strongly  on  the 
SGS  scalar.  The  SGS  scalar  is  also  characterized  by  the  existence  of  diffusion-layer-like  structure 
(ramp-cliffs).  These  characteristics  are  similar  to  the  scalar  PDF  in  the  early  stages  of  initially 
binary  mixing[12].  Therefore,  when  used  as  a  conditioning  variable  for  the  conditional  scalar  FDF, 
the  SGS  scalar  variance  can  reveal  the  structure  of  the  SGS  scalar  that  is  otherwise  averaged  out, 
providing  important  new  SGS  physics.  This  property  of  the  SGS  variance  is  partly  because  it  is  an 
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Figure  1:  Conditional  FDF  of  mixture  fraction,  Z ,  in  the  jet.  The  FDF  is  Gaussian  and  bimodal  for 
small  and  large  SGS  variance  values,  respectively.  The  beta  model  is  able  to  predict  the  Gaussian 
FDF  but  not  the  bimodal  FDF.  The  SGS  variance  normalized  by  its  mean  is  given  in  the  legend. 


Figure  2:  Mean  FJDF  conditional  on  the  SGS  energy  and  the  resolvable-scale  velocity:  (a) 
k'L/(kL)  =  0.33,  (ui)L  =  (u)  (3  m/s),  {u^)l  —  0.  The  FJDF  is  close  to  joint-Gaussian;  (b) 
kL/(kL)  =  4.2,  (u)l  =  (u)  (3  m/s),  (v)i  =  0.  The  FJDF  has  an  approximately  uniform  region. 


important  variable  in  the  inertial-range  dynamics  and  for  characterizing  the  state  of  SGS  mixing. 
Since  the  conditional  scalar  FDF  for  different  SGS  variance  is  dominated  by  different  structures  and 
dynamics,  it  can  potentially  be  modeled  more  accurately  than  the  unconditioned  FDF  (or  PDF), 
leading  to  improved  LES  statistics.  Therefore,  such  conditional  FDFs  are  important  for  studying 
SGS  scalars  and  modeling  SGS  mixing. 

The  mean  FJDF  conditional  on  the  SGS  energy,  and  the  resolvable-scale  velocity,  U£  is 
shown  in  Fig.  2.  For  small  SGS  energy  the  conditional  FJDF  is  close  to  joint-Gaussian,  similar  to 
the  JPDF  in  a  fully  developed  turbulent  jet,  which  is  in  quasi-equilibrium.  For  large  SGS  energy 
{ki/(ki)  =  4.2)  the  FJDF  has  an  approximately  uniform  region,  beyond  which  the  FJDF  decreases 
rapidly  (faster  than  joint-Gaussian).  The  existence  of  the  uniform  region  suggests  that  the  SGS 


Figure  3:  Streamline  plots  of  the  conditional  filtered  viscous  acceleration  for  A/rj  =  125, 
1  ^2)l|(u)l,  kL)  •  (“l )l  =  («i)  and  (u2)l  =  0.  (a)  kL/(kL)  =  1.2, 

(b)  fcf,/ (%)  =  4.6.  The  magnitudes  of  these  variables  are  given  as  grayscales. 


velocity  field  under  this  condition  contains  structures  in  which  the  velocity  varies  approximately 
linearly  in  physical  space.  Examples  of  such  structures  include  axisymmetric  contraction  and 
expansion,  plane  strain,  and  plane  strain  plus  plane  shear.  Since  these  local  structures  occur  at 
large  %,  the  strain  rate  imposed  by  them  is  larger  than  the  average  turbulent  strain  rate  near  the 
filter  scale.  As  a  result,  the  production  of  the  SGS  energy  increases,  and,  when  the  strain  rate  is 
sufficiently  large,  the  “background”  SGS  turbulence  (the  SGS  turbulence  minus  the  structure)  may 
be  undergoing  local  rapid  distortion,  suggesting  that  the  SGS  velocity  field  is  in  non-equilibrium, 
because  rapid  distortion  is  a  linear  process  and  spectral  transfer  is  not  effective  in  responding  to  the 
increased  SGS  energy  production.  Therefore,  the  dissipation  lags  behind  the  production,  resulting 
in  local  non-equilibrium.  Also  see  discussions  in  Part  C. 

Since  wavenumber  information  is  needed  to  describe  turbulence  under  rapid  distortion,  a  con¬ 
ventional  one-point  PDF  model  is  not  sufficient.  To  overcome  this  difficulty,  Van  Slooten  and  Pope 
(1997)  developed  a  velocity-wave-vector  PDF  model  that  gives  the  exact  Reynolds  stress  in  the 
rapid  distortion  limit.  The  present  study  suggests  that  local  rapid  distortion  might  exist  even  when 
the  mean  strain  rate  is  not  large.  Thus  it  may  be  beneficial  to  adopt  this  model  for  LES  so  that 
the  local  rapid  distortion  effects  can  be  taken  into  account. 


In  the  present  study  two  (partial)  components  of  the  viscous  acceleration  were  measured.  The 
conditionally  filtered  viscous  acceleration,  conditional  on  the  velocity 

components  represents  the  transport  velocity  of  the  VFJDF  by  viscous  acceleration  in  the  velocity 
space  and  is  shown  in  Fig.  3  as  streamline  plots.  A  factor  of  1/2  was  used  for  the  u2  component 
because  for  isotropic  turbulence  its  rms  is  twice  that  of  the  u\  component.  The  magnitude  of  the 
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Figure  4:  The  conditionally  filtered  energy  dissipation  rate  ((£|ui,«2)l|(u)l,  ^l)/{^l\^l)  The  con¬ 
ditions  are  the  same  as  in  Fig.  3. 


vector  field,  normalized  by  ( 2kL)1/2/(eL\kL ),  is  given  as  gray  scale  isocontours.  The  streamlines 
generally  flow  towards  a  stagnation  point,  which  is  approximately  at  the  center  of  the  VFJDF.  The 
magnitude  of  the  acceleration  increases  with  the  magnitude  of  the  velocity.  This  dependence  is 
similar  to  that  of  the  scalar  diffusion  on  the  value  of  the  scalar [13,  14,  15].  The  conditional  viscous 

acceleration  in  the  velocity  JPDF  equation,  (i/(-  j-,  2~d^^Ul'U2^'  a^so  ^ias  similar  properties 

[Fig.  11(c)].  The  normalized  magnitude  of  the  conditionally  filtered  viscous  acceleration  increases 
faster  for  larger  the  SGS  energy,  probably  because  the  local  rapid  distortion  at  large  SGS  energy 
causes  the  SGS  fluctuations  and  accelerations  to  increase,  although  the  linear  portion  of  the  SGS 
velocity  does  not  directly  contribute  to  viscous  acceleration.  The  observed  rate  of  increase  might 
also  be  partially  due  to  the  decrease  in  the  normalized  sample-space  variable  v/(2 We  note 
that,  although  the  conditional  viscous  acceleration  is  non-negligible,  the  rms  viscous  acceleration 
becomes  smaller  than  the  rms  pressure  gradient  as  Reynolds  number  increases  [16]. 

The  conditionally  filtered  dissipation  tensor,  {ejj(«i,U2)L>  provides  an  alternative  to  the  condi- 

dv-i  duj 

tionally  filtered  viscous  acceleration  to  close  the  VFJDF  equation,  where  ^  igure  4 

shows  the  surface  plots  for  the  conditionally  filtered  energy  dissipation  rate  (e|ui,«2)z,  normalized 
by  the  filtered  energy  dissipation  rate,  {ei\ fct,}.  In  general,  the  dissipation  increases  with  u\.  For 
small  SGS  energy  the  surface  is  concave,  similar  to  the  dependence  of  the  scalar  dissipation  on 
the  scalar  values  for  small  SGS  variance.  The  surface  plot  for  the  conditional  energy  dissipation 
(e|«i,  U2),  which  appears  in  the  velocity  JPDF  equation,  also  has  a  similar  shape.  For  large  SGS 
energy  there  is  a  slight  bulge  in  the  center  portion  of  the  surface,  suggesting  that  the  local  velocity 
field  responsible  for  the  large  strain  rate  deviates  somewhat  from  a  linear  field  because  the  dissi¬ 
pation  is  constant  in  such  a  field.  Nonetheless,  the  deviation  is  small,  as  indicated  by  the  uniform 
portion  of  the  conditional  VFJDF.  The  moderate  dependence  of  the  dissipation  on  the  SGS  velocity 
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Figure  5:  Conditional  production  of  SGS  energy  (energy  transfer  rate)  normalized  by  the  filtered 
energy  dissipation  rate,  ((( UiUj)L  —  (ui)  L(uj)  L)sij\ki)  /  {£L\kL,)  ■  The  production  exceeds  the  dissi¬ 
pation  for  large  values  of  &£,  indicating  non-equilibrium  SGS  velocity.  The  filter  sizes  are  given  in 
the  legend. 

suggests  that  the  current  practice  of  modeling  the  dissipation  independent  of  the  velocity  is  largely 
justified. 

C.  Structure  of  the  SGS  velocity 

The  non-universal  VFJDF  raises  questions  about  the  structure  of  the  SGS  velocity  field.  We 
investigated  the  structure,  which  is  important  for  modeling  VFJDF.  In  turbulence  under  rapid 
distortion,  the  spectral  energy  transfer  rate  is  small  compared  to  the  production  due  to  the  mean 
strain.  Thus,  the  turbulent  fluctuations  gain  more  energy  than  the  spectral  transfer  (and  than  the 
dissipation)  can  remove,  suggesting  non-equilibrium  turbulence.  The  SGS  turbulence  under  rapid 
distortion  is  also  found  to  be  in  non-equilibrium[17].  To  determine  whether  the  SGS  turbulence 
under  local  rapid  distortion  is  in  non-equilibrium,  we  compared  the  production  rate  of  the  SGS 
energy,  TijSij,  to  the  filtered  dissipation  rate  e^.  Here  is  defined  as  ( U{Uj ) l  —  ( Ui)i,{uj ) l-  We  used 
IfCniSil  +2tx2Si2  +  'T22s22)  and  5^[(|^-)2  +  (f^)2]  as  surrogates  for  the  full  energy  production  rate 
and  the  full  dissipation  rate,  respectively.  Figure  5  gives  the  ratio  of  the  conditional  production 
of  the  SGS  energy  to  the  conditional  dissipation  conditional  on  the  SGS  energy.  Indeed,  the 
conditional  production  exceeds  the  conditional  dissipation  when  the  SGS  energy  is  larger  than  its 
mean  value,  indicating  that  the  spectral  transfer  is  not  effective  in  responding  to  the  increased  SGS 
energy  production  and  that  the  SGS  turbulence  is  in  non-equilibrium.  This  result  is  consistent 
with  those  for  SGS  turbulence  under  rapid  distortion  by  the  mean  flow[17].  The  non-equilibrium 
of  the  SGS  velocity  is  similar  to  the  non-equilibrium  of  the  SGS  scalar  when  the  SGS  variance  is 
large[5].  Like  the  non-equilibrium  SGS  scalar,  which  contains  scalar  structures  (diffusion  layers), 
the  non-equilibrium  SGS  velocity  also  contains  structures. 
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Figure  6:  The  parameters  B  as  a  function  of  ki/ei  for  A/77  =  125  for  different  ki  values.  The 
collapse  of  the  curves  indicates  that  the  VFJDF  can  largely  be  characterized  by  k^ / / {{ki) / (ei)) . 


To  further  examine  the  effects  of  the  non-equilibrium  SGS  turbulence  on  the  VFJDF  we  com¬ 
puted  the  conditional  VFJDF  with  the  filtered  dissipation  ei,  used  as  an  additional  conditioning 
variable.  The  SGS  dissipation  time  scale,  fc^/ejr ,  can  be  used  as  an  indicator  of  the  degree  of  the 
non-equilibrium:  large  SGS  time  scale  indicates  that  dissipation  is  not  fast  enough,  and  the  SGS 
turbulence  is  in  non-equilibrium. 

To  quantify  the  dependence  of  the  shape  of  the  VFJDF  on  ki  and  e/,,  we  used  two  slices  of 
the  VFJDF  going  through  the  origin  in  the  v\  and  vi  directions,  respectively.  We  computed  the 
conditional  kurtosis  and  skewness  of  the  VFJDF,  KU1 \U2  and  SU1\U2  (for  =  0),  which  can  be 
obtained  from  the  conditional  fourth  and  third  moments  of  {ful\U2(vl)\^L^L,  {u)l),  normalized 
by  the  second  moment  raised  to  the  proper  power,  where  fm\U2  is  the  FDF  of  the  u\  component 
conditional  on  the  u 2  component  (essentially  a  slice  in  the  u\  direction).  The  conditional  kurtosis 
and  skewness  KU2\U1  and  SU2|U1  can  also  be  obtained  in  the  same  way.  We  combined  K  and  S 
to  form  a  parameter  B  =  K  —  S2,  proposed  by  Atkinson[18].  This  parameter  takes  the  value 
of  unity  for  a  double-delta  PDF  regardless  of  its  symmetry,  therefore,  is  a  better  measure  of  the 
intermittency  (or  lack  of)  than  the  kurtosis.  For  a  uniform  distribution  B  has  a  value  of  1.8.  It 
has  been  used  to  characterize  the  bimodal  scalar  FDF [4]  and  velocity  increment  PDF[19]. 

We  examined  BU1  \U2  for  u"  =  0  and  BU2 |Ul  for  u'[  =  0.  If  the  VFJDF  is  perfectly  uniform  both 
parameters  have  the  value  of  1.8.  Figure  6  shows  that  B  generally  decreases  with  ei,  consistent 
with  the  observation  that  the  VFJDF  changes  towards  a  uniform  shape  as  decreases.  Due  to 
the  opposite  trends  of  B  with  respect  to  fcj,  and  e^,  we  plotted  BU1  \U2  and  BU2  |U1  against  ki/ei 
for  three  ki  values.  Figure  6  shows  that  both  BU1\U2  and  BU2 jui  decrease  as  ki/e^  increases. 
Furthermore,  the  curves  for  BU2\U1  appear  to  collapse  for  different  k^  values,  although  the  results 
for  BU1  \U2  still  show  some  kj_,  dependence.  Nonetheless,  these  results  indicate  that  k^/t^  alone 


13 


Figure  7:  Dependence  of  the  bimodality  parameter  on  k?  jtT .  The  averaging  domain  sizes  axe  given 
in  the  legend.  The  values  of  p  is  approximately  0.76.  The  the  dashed  and  solid  lines  represent 
kT/(kr)  values  of  1.0  and  4.0  respectively.  Except  for  r/p  =  40,  B&u  is  essentially  determined  by 
k.P/er  and  is  independent  of  kr. 

can  largely  characterize  the  shape  of  the  VFJDF.  Specifically,  the  VFJDF  is  uniform  as  long  as 
ki/f-L  is  large,  i.e.,  when  the  SGS  velocity  is  in  strong  non-equilibrium,  regardless  of  the  value 
of  ki,.  This  condition  is  less  restrictive  than  that  of  &£,/(£/,)  S>  1.  The  observed  importance  of 
kh/f-L  in  determining  BU1\U2  and  BU2\U1  suggests  that  the  degree  of  non-equilibrium  is  important 
for  determining  the  distribution  of  the  SGS  velocity  field. 

To  further  understand  the  structure  of  the  SGS  velocity,  we  also  examined  the  conditional 
velocity  increment  with  separations  the  same  as  the  filter  sizes.  The  increment  explicitly  contains 
(two-point)  spatial  information  whereas  the  VFJDF  only  contains  information  implicitly.  Therefore, 
conditional  velocity  increment  statistics  provide  essential  information  complementing  the  VFJDF. 

In  Fig.  7  we  plotted  the  bimodality  parameter  conditional  on  k l  and  €i  (here  kr  and  er  are 
used  for  ki  and  ex,  to  indicate  that  the  results  are  for  increments)  and  for  the  velocity  increment 
PDF  as  a  function  of  fc^/cjr,,  where  p  is  the  scaling  exponent  for  (eT\kr)  ~  kf.  We  choose  the 
ratio  of  (er\kr)/er  ~  fcf/er  as  the  variable  because  it  is  an  indicator  of  the  deviation  of  er  from  its 
conditional  mean.  For  a  randomly  oriented  plane  strain  field,  B&u  is  bimodal.  For  all  separations 
considered,  B&u  decreases  with  fcf?/er,  and  the  increment  PDF  is  bimodal  for  large  kP/er  values. 
Furthermore,  B&u  is  approximately  independent  of  the  kr  values,  indicating  that  the  bimodality 
is  largely  determined  by  the  relative  magnitude  of  er  to  its  conditional  mean,  suggesting  that  the 
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bimodality  is  closely  related  to  the  degree  of  non-equilibrium  of  the  local  velocity  field  because  as 
kp/er  increases  the  production  of  kr  increases  relative  to  the  dissipation.  The  bimodality  parameter 
Bav  has  a  similar  trend,  but  there  is  still  some  weak  dependence  on  kr.  Furthermore,  it  appears 
that  the  curves  for  different  filter  scales  can  be  collapsed  by  shifting  them  horizontally,  i.e.,  by  a 
multiplying  factor  to  kp/eT.  Except  for  r/77  =  40,  this  factor  ranges  from  1.1  to  1.5.  Thus,  the 
need  for  this  factor  might  be  a  result  of  the  limited  Reynolds  number  (and  the  scale  separations 
between  r  and  the  integral  length  scales)  in  the  present  study.  At  higher  Reynolds  numbers,  we 
might  expect  improved  scaling  regions  for  k~p(er\kr)  and  collapse  of  the  results  for  the  bimodality 
parameter.  This  result  suggests  that  there  may  be  “universal”  dependence  of  the  bimodality  on 
kp/er.  Thus,  the  results  provide  further  evidence  of  the  plane  strain  structure  inferred  from  the 
VFJDF. 

D.  Velocity-scalar  FJDF,  conditionally  filtered  scalar  diffusion,  and  dissipation 

This  part  focuses  on  the  coupling  between  the  SGS  velocity  and  scalars,  which  is  essential  for 
modeling  the  effects  of  the  velocity  on  the  SGS  mixing. 

The  mean  FJDF  conditional  on  four  variables:  the  SGS  energy,  the  resolvable-scale  velocity, 
the  SGS  scalar  variance,  and  the  resolvable-scale  scalar  has  qualitatively  different  forms  depending 
on  the  conditioning  variables.  For  small  kjj  (<  (ki))  and  {(fi"2) l  (<  (0"2)),  the  conditional  FJDF 
is  close  to  joint-Gaussian  (not  shown),  similar  to  the  JPDF  in  a  fully  developed  turbulent  jet, 
which  is  generally  in  quasi-equilibrium.  For  large  ( <$>"2)l  the  FJDF  is  bimodal  at  4>"/{4>,,2)l  ~  ±1 
regardless  of  the  values  of  fcjr,  (Fig.  8),  consistent  with  a  scalar  diffusion  layer  structure.  For  small  ki 
[Fig.  8(a)]  the  correlation  between  the  SGS  scalar  and  the  velocity  component  u\  is  low,  indicating 
small  conditional  SGS  flux.  For  large  ki  [Fig.  8(b)]  the  correlation  is  strong,  and  the  two  peaks 
of  the  FJDF  are  also  separated  in  the  velocity  space  with  a  “neck”  near  u\  =  0.  Thus,  scalar 
mixing  is  expected  to  be  most  intense  near  this  point  and  depends  strongly  on  the  velocity,  further 
demonstrating  the  need  for  mixing  models  to  depend  on  the  velocity.  This  FJDF  shape  is  consistent 
with  a  scalar  diffusion  layer  associated  with  a  local  plane  strain  (converging-diverging)  velocity  field. 
In  such  a  structure  the  scalar  gradient  is  largest  in  the  diffusion  layer,  where  the  SGS  velocity  is 
the  smallest.  The  structure  is  similar  to  the  counter-flow  model  for  laminar  flamelets.  The  FJDF 
peaks  in  Fig.  8(a),  on  the  other  hand,  overlap  in  the  velocity  space.  Thus  mixing  can  occur  over  a 
wide  range  of  SGS  velocity  and  is  less  dependent  on  the  velocity. 

In  the  present  study  the  scalar  diffusion  and  dissipation  were  obtained  using  the  streamwise 
derivatives.  For  small  {4>"2)l  the  diffusion  generally  has  linear  dependence  on  both  the  SGS  velocity 
and  scalar.  Figure  9(a)  shows  the  conditionally  filtered  diffusion  for  small  ki,  (0.52 (ki)).  In  the 
center  portion  the  isocontours  can  be  approximated  by  equally-spaced  straight  lines,  consistent 
with  a  previous  linear  model.  For  large  ki  the  isocontours  have  a  larger  slope,  indicating  a  stronger 
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Figure  8:  Mean  FJDF  conditional  on  ki,  (u )l,  (ft'2)L,  and  (4>)l  on  the  jet  centerline  for  («i)z,  = 
<«!)  (3.07  m/s)  and  (uj)L  =  0:  (a)  kL/{kL)  =  3.04  and  (ft'2)L/(ft'2)  =  6.09;  (b)  kL/(kL)  =  0.64 
and  (ft'2)L/(ft'2)  —  6.09.  The  bimodal  FJDF  in  (b)  is  strongly  asymmetric  in  u\,  indicating  large 
SGS  flux  and  suggesting  that  mixing  occurs  primarily  near  u\  =  0. 


Figure  9:  Mean  conditionally  filtered  scalar  diffusion:  (a)  kiKkft)  =  0.8  and  (ft'2)L/(ft'2)  =  0.94; 
(b)  kjJ(ki)  =  0.52  and  (ft'2)  l  I  (ft'2)  =  6.09.  The  diffusion  in  (a)  is  can  be  approximated  using  a 
previous  linear  model  (straight  lines)  whereas  that  in  (b)  has  a  S-shaped  region. 
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Figure  10:  Normalized  mean  conditionally  filtered  scalar  dissipation:  (a)  kjJ{ki)  —  0.64  and 
(f2)L/(r2)  —  4.67;  (b)  ki/{ki)  =  2.43  and  (<£//2)z,/(<£//2)  =  4.67.  The  dissipation  in  (b)  shows 
that  mixing  occurs  primarily  near  u\  =  0. 

dependence  on  u\.  The  results  are  consistent  with  velocity-scalar  fields  having  a  joint  normal  JPDF. 

For  large  SGS  variance  the  surface  plot  of  the  diffusion  is  S-shaped  and  depends  on  uj  for 
certain  ranges  of  &£,  values.  Figure  9(b)  shows  the  conditionally  filtered  diffusion  for  large  (4>"2)l 
(6.09(^"2))  but  small  k^.  In  addition  to  an  approximately  linear  trend,  an  S-shaped  dependence 
on  (f>  exists  for  all  u\  values,  consistent  with  the  FJDF  in  Fig.  8(a),  which  shows  that  the  bimodal 
scalar  distribution  spans  a  wide  range  of  velocity.  For  large  fcjr,  the  S-shaped  region  appears  to  be 
limited  to  near  iq  =  0  (not  shown),  suggesting  that  large  magnitudes  of  diffusion  occur  at  small 
SGS  velocity.  This  issue  is  further  discussed  in  the  following,  along  with  the  dissipation. 

The  conditionally  filtered  scalar  dissipation  provides  an  alternative  to  the  conditionally  filtered 
scalar  diffusion  to  close  the  FJDF  equation.  For  small  the  dissipation  generally  depends 

weakly  on  both  the  velocity  and  scalar  (not  shown),  consistent  with  the  approximately  joint  nor¬ 
mal  FJDF.  For  small  ki  and  large  {4>',2)l  the  dissipation  shown  in  Fig.  10(a)  has  a  bell-shaped 
dependence  on  <f>  but  is  weakly  dependent  on  iq.  The  dependence  on  <j)  is  consistent  with  the  FJDF 
in  Fig.  8(a)  and  the  S-shaped  diffusion  plot  in  Fig.  9(b).  For  large  k^  and  (<f)"2) l  the  dissipation 
is  largest  at  <f>  =  0  and  u\  =  0  [Fig.  10(b)]  and  decreases  as  velocity  increases,  indicating  a  strong 
dependence  on  the  SGS  velocity,  consistent  with  the  FJDF  in  Fig.  8(b)  and  the  notion  that  the 
diffusion  layer  is  associated  with  a  converging-diverging  (e.g.,  plane  strain)  velocity  field.  In  such  a 
structure  the  velocity  component  normal  to  the  diffusion  layer  is  small  near  the  layer.  This  result 
also  suggests  that  the  S-shaped  surface  for  the  scalar  diffusion  is  limited  to  small  u\  values.  Since 
diffusion  layers  have  been  shown  to  exist  at  very  high  Reynolds  numbers,  the  observed  dependence 
of  scalar  diffusion  and  dissipation  on  the  velocity  is  expected  to  be  Reynolds-number  independent. 

This  study  showed  that  there  is  a  strong  statistical  dependence  between  scalar  mixing  (dissipation- 
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scale  scalar)  and  the  SGS  velocity  field.  The  scalar  and  velocity  structure  for  large  SGS  variance 
(especially  when  the  SGS  energy  is  also  large)  is  similar  to  the  structure  of  laminar  flamelets.  Thus 
the  ability  of  mixing  model  to  reflect  such  a  structure  and  the  localness  of  mixing  is  important. 
The  scalar  diffusion  conditional  on  the  SGS  velocity  under  such  conditions  is  expected  to  be  non¬ 
zero  and  independent  of  Reynolds  number.  The  results  further  demonstrate  the  importance  of  the 
velocity  in  mixing  and  of  including  velocity  in  mixing  models. 

E.  The  joint  structure  of  SGS  velocity  and  scalar 

The  velocity-scalar  FJDF  observed  suggests  coupled  structures  of  the  SGS  velocity  and  scalars. 
We  investigated  the  velocity-scalar  increment  JPDFs  to  further  examine  the  structures. 

When  conditioned  on  er  and  Xr  (again,  r  is  used  in  the  subscript  instead  of  a  L  to  indicate  that 
the  results  are  for  increments),  the  joint  PDF  of  velocity  and  scalar  increments,  PAv,A<t>\er,xr  c^ose 
to  joint  Gaussian  for  a  range  of  er  and  \r  values,  consistent  with  Kolmogorov’s  refined  similarity 
hypotheses  (K62)  [20].  When  conditioned  on  the  local  kinetic  energy  and  the  local  scalar  variance, 
PAv,A4>\^rM"2)r  ^ias  different  shapes  depending  on  the  values  of  the  conditioning  variables.  For  small 
(<f>"2)r  and  Ay,  it  is  close  to  joint  normal  [Fig.  11  (a)].  For  small  (<^"2)r  but  large  Ay,  it  is  bimodal  in 
Av ,  indicating  a  well-mixed  local  scalar  field  but  a  nonequilibrium  local  velocity  field.  The  results 
in  Part  C.[21]  have  shown  that  when  Ay  is  large  the  local  turbulence  is  under  rapid  distortion, 
and  the  local  velocity  field  contains  plain  strain.  Figure  11  (b)  suggests  that  the  plane  strain  does 
not  alter  the  distribution  of  a  well-mixed  scalar.  For  large  ( <j>"2)r  but  small  Ay,  P/^v,A<t>\kr,(<t>"2)r 
is  bimodal  in  A(f>,  indicating  a  nonequilibrium  local  scalar  field  but  an  equilibrium  local  velocity 
field.  Although  the  velocity  is  in  quasi-equilibrium  at  scale  r,  the  bimodal  JPDF  suggests  that 
there  may  be  nonequilibrium  velocity  fields  at  smaller  scales.  Therefore,  the  diffusion  layers  under 
such  conditions  might  be  highly  wrinkled.  When  both  ( <t>"2)r  and  Ay  are  large,  PAv,&<j>\kr,(<l>"2)r  *s 
quad-modal,  with  one  peak  in  each  quadrant  of  the  Av  —  A <f>  plane.  This  JPDF  shape  has  not  been 
previously  observed  for  inertial-range  turbulence  and  provides  further  evidence  of  the  combination 
of  a  local  plane  strain  and  ramp-cliff  (diffusion  layer)  structure. 


The  similarity  between  the  SGS  scalar  at  large  SGS  variance  and  binary  mixing  suggests  that  the 
SGS  scalar  is  in  nonequilibrium.  To  further  examine  the  role  played  by  spectral  nonequilibrium,  we 

rv/  j  v 

computed  the  conditional  production  rate  of  ( <f>"2)r ,  —  ~{{uj4’)r  —  (uj)r(<P)r)—^ — conditional 

on  (4>"2)r  and  Xr,  respectively  (Fig.  12).  The  averaging  domain  size  r  is  125 rj.  When  conditioned 
on  (4>''2)r,  the  conditional  mean  of  V^T  normalized  by  (Xr|(<^,/2)r)  is  small  for  {4>"2)r/ <  1 
and  increases  with  (<p"2)T  (~  (<^"2)°'8),  indicating  that  the  local  scalar  field  changes  from  quasi¬ 
equilibrium  to  nonequilibrium  as  {4>n2)T  increases.  (A  decaying  scalar  field  has  similar  characteristics 
to  a  quasi-equilibrium  scalar  field).  On  the  other  hand,  when  conditioned  on  the  scalar  dissipation, 
P<f>r/Xr  only  increases  by  a  factor  of  two  over  two  orders  of  magnitude  increase  in  Xr ■  The  results 
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Figure  11:  Conditional  JPDF  of  Av  i 
and  (< f2)r /(<!>' /2>  are  (a)  0.06,  0.3;  (b) 


Figure  12:  Conditional  mean  of  the  local  scalar  variance  production  conditional  on  (4>"2)r  and 
Xr,  respectively.  The  averaging  domain  size  r  is  125r?. 


Figure  13:  Dependence  of  the  bimodality  parameter  on  {<t>"2)T/Xr-  The  values  of  ((f>"2)r  /  (<t>'/2)  are 
given  in  the  legend.  B  is  essentially  determined  by  (</>"2)r/Xr,  independent  of  (<f>"2)T- 


suggest  that  when  conditioned  on  Xr,  the  conditional  inertial-range  scalar  is,  on  average,  in  quasi¬ 
equilibrium,  whereas,  when  conditioned  on  { <j>"2)r ,  the  degree  of  nonequilibrium  of  the  conditional 
inertial-range  scalar  increases  with  (< p"2)r ■  These  results  provide  further  evidence  that  the  close- 
to-Gaussian  conditional  increment  PDFs  are  a  result  of  quasi-equilibrium  scalar  fields  as  implied 
in  Kolmogorov’s  refined  hypothesis,  and  suggests  that  the  bimodal  conditional  FDFs  axe  closely 
related  to  the  nonequilibrium  of  the  conditional  inertial-range  scalar. 

To  quantify  the  dependence  of  the  conditional  PDFs  on  the  degree  of  nonequilibrium,  we  com¬ 
puted  the  bimodality  parameter [18],  B  =  K  —  S2,  as  a  function  of  {4>"2)T/Xr  for  three  {<j>"2)T  values. 
Figure  13  shows  that  B  is  approximately  independent  of  {<t>"2)r  and  is  largely  a  function  of  (4>"2)r/Xr 
alone.  Therefore,  the  conditional  PDF  is  largely  determined  by  the  degree  of  nonequilibrium.  For 
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Figure  14:  Conditional  mean  of  the  locally  averaged  energy  dissipation  (er|A:r)  multiplied  by  k~p, 
where  p  is  the  scaling  exponent  of  (er  |fcr).  Approximate  scaling  ranges  exist  for  all  averaging  domain 
sizes  (given  in  the  legend). 

large  (4>"2)r/Xr  values  B  appears  to  decrease  monotonically  towards  unity  (the  slight  increase  of  B 
for  very  large  ( <f>"2)r/Xr  is  due  to  insufficient  statistical  convergence  there),  indicating  that  the  con¬ 
ditional  increment  PDF  approaches  a  double  delta  function  for  strong  nonequilibrium  local  scalar 
fields.  The  results  suggest  that  the  degree  of  spectral  nonequilibrium  can  be  parameterized  by  the 
ratio  (</>"2)r/Xr  for  three  (4>"2)r,  which  can  be  modeled  in  LES. 

F.  Issues  of  modeling  the  filtered  energy  and  scalar  dissipation  rates 

In  the  velocity-scalar  FJDF  approach,  the  filtered  energy  dissipation  and  the  filtered  scalar 
dissipation  require  modeling.  The  current  practice  uses  the  spectral  equilibrium  assumption  and 
therefore  are  inconsistent  with  the  results  discussed  above.  We  investigated  the  issues  of  their 
modeling. 

We  first  examined  the  statistical  relationship  between  kr  and  er.  Using  Kolmogorov’s  refined 
hypotheses[20],  one  can  write  kr  =  V^er^r2/3  for  an  inertial-range  separation  r,  where  Vi  is  a 
universal  stochastic  variable  independent  of  er.  Thus,  the  conditional  mean  (kr\er)  scales  as  e^3, 
which  has  been  confirmed  experimentally [22].  However,  the  scaling  exponent  for  the  reversed 
conditional  mean,  (er|fcr)  ~  kp,  does  not  follow  the  K62  prediction  of  p  =  3/2[22]  and  thus  is 
“anomalous”.  Here  we  obtained  the  scaling  exponents  p  by  multiplying  k~p  to  (er|fcr)  to  obtain  a 
flat  portion  of  the  curve  (Fig.  14).  For  r/rj  =  40,  there  is  a  clear  scaling  range  between  kr/(kr)  = 
0.01  and  3  with  a  scaling  exponent  of  0.763.  For  larger  separations  there  are  approximate  scaling 
ranges  where  the  variations  of  k~p(er\kr)  are  within  30%  of  the  peak  values.  The  exponents  obtained 
are  between  0.765  and  0.77.  We  might  expect  that  at  higher  Reynolds  numbers  more  definitive 
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scaling  regions  will  emerge  for  these  separations  as  well. 

A  possible  explanation  for  the  deviations  of  the  observed  scaling  exponent  p  from  predictions 
of  Kolmogorov’s  refined  hypotheses  is  that  the  hypotheses  are  based  on  quasi-equilibrium  inertial- 
range  spectral  energy  transfer,  that  is,  the  dissipation  balances  (or  is  larger  than)  the  spectral 
transfer.  Under  such  conditions  the  local  energy  dissipation  time  kr/er  scales  as  the  local  energy 
production  (spectral  transfer)  time  r/fc^2.  When  er  is  used  as  the  conditioning  variable,  this  quasi¬ 
equilibrium  condition  is  likely  to  be  satisfied  because,  for  a  given  magnitude  of  energy  dissipation, 
there  generally  is  a  corresponding  amount  of  the  spectral  transfer.  Thus  (kr  |er)  scales  as  e2^3r2/3, 
but  the  reverse  scaling,  (er\ kr)  ~  /c^2,  is  not  necessarily  true  because  an  increase  in  the  spectral 
transfer  near  the  scale  r  can  result  in  large  kr  and  a  non-equilibrium  local  velocity  field  with  the 
dissipation  lagging  behind.  The  nonequilibrium  spectral  transfer  results  in  a  slower  increase  of  er 
than  the  equilibrium  prediction  of  er  ~  .  In  fact,  the  observed  scaling  er  ~  k?  shows  that 

the  ratio  of  the  local  energy  dissipation  time  scale  to  the  local  production  time  scale  varies  as 
£3/2-p/r,  which  increases  with  kr.  Thus,  it  appears  that  the  anomalous  scaling  is  due  to  the  non¬ 
equilibrium  of  the  local  velocity  fields,  The  results  show  that  models  for  e/,  must  take  into  account 
the  nonequilibrium  spectral  transfer. 

Next,  we  examined  the  issue  of  modeling  XL ■  To  understand  the  variations  of  the  degree  of 
non-equilibrium  of  the  local  scalar  fields,  it  is  useful  to  examine  the  relationship  between  (<f>"2)r  and 
Xr-  In  Fig.  15  we  plot  two  conditional  means,  ((<£"2)r|Xr)  and  (Xr\(<f>"2)r)-  The  former  increases  as 
Xr  (the  results  are  similar  when  Xr^r  ^3  is  used  as  the  conditioning  variable),  consistent  with  the 
K62  prediction 

(4>"\  =  ^2Xr€r"1/3r2/3,  (13) 

where  is  also  a  universal  stochastic  variable  independent  of  both  Xr  and  er.  However,  the  latter 
(as  well  as  (xr^r  1^3|(0//2)r),  not  shown)  scales  approximately  as  (<^>,/2)°-7;  therefore  it  is  at  variance 
with  the  K62  prediction  of  ~  ((f>"2) r.  A  possible  reason  for  the  deviation  of  (Xr\(<t>"2)r)  from  the 
K62  scaling  is  that  the  conditional  inertial-range  scalar  is  in  quasi-equilibrium  when  conditioned 
on  Xr,  whereas,  when  conditioned  on  (<f>"2)r  it  can  be  in  nonequilibrium:  For  a  given  value  of  Xr 
there  generally  is  a  corresponding  amount  of  spectral  transfer  (scales  as  (^//2)r/r^r).  Therefore,  the 
scalar  is  in  quasi-equilibrium.  On  the  other  hand,  when  (<f>"2)r  (therefore  the  spectral  transfer)  is 
increased,  Xr  generally  lags  behind  since  it  takes  time  for  the  increased  spectral  transfer  rate  to  reach 
the  dissipation  scales,  resulting  in  non-equilibrium  local  scalar  fields  and  a  smaller  scaling  exponent 
f°r  (Xr|(</>,,2)r)-  Therefore,  models  for  xl  must  contain  the  effects  of  spectral  nonequilibrium. 

The  influence  of  the  velocity  field  is  also  important  for  modeling  XL  and  ( <f>"2)L  (the  latter  is 
needed  in  non-PDF-based  LES)  and  was  examined.  The  conditional  mean  of  (4>"2)l  conditional  on 
the  filtered  energy  dissipation  ei  =  (e)z,  and  the  SGS  energy,  ((4>"2)lVl,  ^l),  is  given  in  Fig.  16(a). 
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Figure  15:  Conditional  mean  of  the  locally  averaged  scalar  dissipation  ((<f>"2)r\Xr)  and  the  reversed 
conditional  mean  {Xr\{<t>"2)r)- 


Figure  16:  Conditional  SGS  scalar  variance  (a),  SGS  variance  production  (b),  and  filtered  dissipa¬ 
tion  rate  (c)  conditional  on  and  ki,. 

The  results  with  sufficient  statistical  convergence  are  limited  to  an  elliptic  region  that  approximately 
corresponds  to  the  isocontours  of  the  e/,  —  k^  JPDF[21].  At  e£/(e)  =  1  and  &z,/(A;jr,)  =  1  the  SGS 
variance  is  close  to  its  mean  value  {{{(j)"2)  l\^Li  &l)  —  0.9(</>"2)),  indicating  that  the  SGS  scalar  is 
close  to  its  average  conditions.  Toward  the  lower  right  portion  of  the  e^  —  ki  plane  the  SGS  variance 
increases. 

To  understand  this  trend  we  computed  the  conditional  means  of  the  SGS  variance  production, 
(P^lI^Li  &£<)•  The  result  given  in  Fig.  16(b)  is  normalized  by  (xil|6L,^l)  and  therefore  repre¬ 
sents  the  balance  between  the  production  and  dissipation  and  the  degree  of  nonequilibrium  of  the 
SGS  scalar.  Similar  to  the  SGS  variance,  the  production  also  increases  toward  the  lower  right 
portion  of  the  ei  —  ki  plane.  In  addition,  the  isocontours  have  approximately  the  same  slopes 
as  those  of  ({<//,2)l|€l,  /c^),  suggesting  that  ({(f)"2)  lVl,  ^l)  is  largely  determined  by  the  degree  of 
non-equilibrium  of  the  SGS  scalar  and  increases  with  it. 
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The  conditionally  filtered  scalar  dissipation  rate,  (xil|cL)^l)>  has  isocontours  approximately 
perpendicular  to  those  of  {(4>"2) [Fig.  16(c)].  Therefore,  xiL  changes  the  fastest  when  the 
degree  of  nonequilibrium  of  the  SGS  scalar  is  held  constant  (along  the  isocontours  of  ((<^,,2)l|cl»  &£,) 
and  &l)).  For  quasi-equilibrium  spectral  transfer,  the  filtered  dissipation  scales  as  XL  ~ 

where  t^l  is  the  SGS  scalar  time  scale.  Therefore,  Fig.  16(c)  suggests  that  the  increase 
of  Xii  along  the  isocontours  of  (4 >"2)z,  is  due  to  the  decrease  in  t^l  caused  by  the  velocity  field. 
The  results  in  Figs.  16(b)  &  (c)  suggest  that  the  SGS  velocity  causes  (<f>"2) l  and  xiL  to  increase  by 
increasing  the  degree  of  non-equilibrium  of  the  SGS  scalar  and  by  decreasing  the  SGS  scalar  time 
scale,  respectively. 

To  further  examine  the  effects  of  the  velocity  field,  we  considered  the  equilibrium  property 
of  the  velocity  field.  A  measure  of  the  degree  of  nonequilibrium  of  the  SGS  velocity,  which  can 
be  obtained  from  the  data,  is  the  conditional  production  of  the  SGS  energy  (Pi\ki,  cl)Al,  where 
Pi  =  —({uiui)i  —  (u\)i(ui)i)dui/dx\,  i  =  1  and  2.  The  isocontours  of  Pi  (not  shown)  in  the  region 
ki  >  ( kL )  have  slopes  similar  to  those  of  £l}/xil-  The  results  suggest  that  the  observed 

trends  for  ((4>''2)i\ei,  ki)  and  (xil|cl,  ki)  are  largely  related  to  the  degree  of  nonequilibrium  and 
the  time  scale  of  the  SGS  velocity.  Therefore,  the  velocity  affects  {{4>"2)i\ei,  ki)  and  (xiL\£L,ki) 
through  different  processes  that  appear  to  have  “orthogonal”  effects. 

These  results  have  implications  for  SGS  modeling.  The  filtered  scalar  dissipation  has  been 
modeled  by  assuming  local  equilibrium  XL  =  P<j>L[ 23].  Figure  16(b)  suggests  that  the  model  can  be 
improved  by  including  ei  and  kL  to  account  for  the  departure  from  equilibrium  spectral  transfer. 
The  production  P^i  is  available  in  velocity-scalar  FJDF  calculations  and  can  be  used  to  model 
Xl-  In  addition,  for  LES  that  solves  scalar  transport  equations,  P^l/xl  can  be  modeled  using 
({(/>"2)L\eL,ki)  and  a  scalar  time  scale  determined  by  ejr,  and  ki .  The  results  for  ((0//2)z,|e£,fcz,) 
can  also  be  used  to  test  models  for  ( 4>n2)i ■ 

G.  Preliminary  study  of  the  SGS  mixing  of  the  mixture  fraction  and  the  SGS  flame 
structure  in  turbulent  flames 

Using  the  knowledge  of  SGS  mixing  gained  in  non-reacting  flows  we  began  to  conduct  a  prelimi¬ 
nary  study  of  the  mixture  fraction  filtered  mass  density  function  (FMDF),  the  SGS  mixture  fraction 
structure,  and  their  influences  on  the  SGS  flame  structure  in  turbulent  flames.  Experimental  data 
obtained  in  Sandia  flames  D  and  E  by  Karpetis  and  Barlow  (www.ca.sandia.gov/TNF/abstract.html) 
were  used  to  analyze  the  SGS  mixture  fraction  and  the  flame  structure.  One-dimensional  filtering 
in  the  radial  direction  on  data  segments  was  employed  to  compute  the  FMDF  and  other  resolvable- 
scale  and  subgrid-scale  variables. 

Conditional  FMDF  conditional  on  the  Favre-filtered  scalar  (mixture  fraction)  and  SGS  scalar 
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'  *  (a)  «  (b) 

Figure  17:  The  mean  conditional  scalar  FMDF  (a)  and  the  mean  conditionally  filtered  scalar 
dissipation  (b)  obtained  in  Sandia  flame  D  at  x/D=15  with  a  filter  size  of  A  =  3.0mm.  The  filtered 
mixture  fraction  (£)l  is  set  to  the  stoichiometric  mixture  faction,  £s,  and  the  SGS  variance  values 
(£,,2)i  are  given  in  the  legend.  Both  the  FMDF  and  the  conditionally  filtered  scalar  dissipation 
indicate  qualitatively  different  SGS  mixing  regimes  and  SGS  mixture  fraction  structure  for  small 
and  large  variance  values. 

variance  shows  a  similar  trend  to  the  non-reacting  results.  For  SGS  scalar  variance  small  compared 
to  its  mean,  the  FMDF  is  not  far  from  Gaussian,  and  the  conditionally  filtered  scalar  dissipation 
rate  depends  weakly  on  the  SGS  scalar,  indicating  that  the  SGS  scalar  is  well  mixed  (Fig.  17). 
For  large  SGS  variance,  however,  the  FMDF  becomes  bimodal,  and  the  conditionally  filtered  scalar 
dissipation  is  bell-shaped,  indicating  the  existence  of  a  diffusion  (dissipation)  layer  structure  that 
is  similar  to  the  mixture  fraction  profile  in  the  counter-flow  model  for  laminar  flamelets.  For  the 
measurement  locations  considered  (up  to  30  jet  diameters  downstream)  the  difference  in  mixture 
fraction  values  for  the  two  FMDF  peaks  is  generally  larger  than  the  reaction  zone  width  in  the  mix¬ 
ture  fraction  space,  and  therefore  the  mixing  field  under  such  conditions  can  support  flamelets.  The 
conditionally  filtered  temperature  near  the  stoichiometric  mixture  fraction  decreases  progressively 
with  increasing  SGS  scalar  variance.  Furthermore,  local  extinction  events  appear  to  occur  mostly 
when  the  SGS  scalar  variance  is  large,  suggesting  the  possibility  of  flamelet  extinction.  These  pre¬ 
liminary  results  suggest  that  the  mixing  regimes  and  the  associated  mixture  fraction  structure  could 
potentially  have  a  strong  impact  on  the  combustion  regime  and  extinction/reignition  in  turbulent 
nonpremixed  flames. 

IV  Conclusions  and  suggestions  for  LES 

The  results  obtained  in  this  research  significantly  advanced  the  understanding  of  the  SGS  ve¬ 
locity  and  scalar,  which  is  important  for  modeling  SGS  mixing  in  LES  of  turbulent  combustion. 
The  issues  include  the  velocity  VFJDF  and  velocity-scalar  FJDF,  the  structure  of  the  SGS  velocity 
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and  scalar,  the  FJDF  transport  equation,  and  the  filtered  energy  and  scalar  dissipation  rates.  The 
specific  conclusions  and  suggestions  for  LES  are: 

•  The  velocity-scalar  FJDF  shows  two  mixing  regimes:  the  equilibrium  spectral  transfer  regime 
with  close-to-Gaussian  velocity  and  scalar  FJDFs,  and  the  nonequilibrium  regime  with  bi- 
modal  scalar  FDF  and  uniform  velocity  FJDF.  The  velocity-scalar  structure  in  the  nonequi¬ 
librium  regime  is  similar  to  that  of  a  laminar  flamelet  model.  Therefore,  consistency  of  mixing 
models  with  the  flamelet  limit  is  important. 

•  The  velocity  FJDF  shows  that  the  SGS  turbulence  is  under  local  rapid  distortion  when 
the  SGS  velocity  is  in  nonequilibrium.  Eddy  viscosity  models  (or  when  an  eddy  viscosity 
is  implied)  are  inconsistent  with  local  rapid  distortion  because  the  evolution  of  turbulence 
under  rapid  distortion  is  determined  by  the  amount  of  strain,  not  the  strain  rate.  The  velocity- 
wavevector  PDF  model  of  Slooten  and  Pope[24]  could  be  adapted  to  LES  to  model  the  local 
rapid  distortion.  However,  this  model  may  significantly  increase  the  computational  cost. 

•  The  velocity-scalar  FJDF  and  the  dependence  of  the  conditionally  filtered  scalar  dissipation 
and  diffusion  on  the  SGS  velocity  further  demonstrate  the  need  to  include  velocity  information 
in  mixing  models  (e.g.,  the  velocity  conditioned  model  of  Fox[25]).  Currently  employed  mixing 
models  are  generally  independent  of  the  velocity. 

•  The  degree  of  spectral  nonequilibrium  can  largely  be  characterized  and  potentially  parame¬ 
terized  by  the  SGS  scalar  variance  and  filtered  scalar  dissipation  for  the  SGS  scalar  and  by  the 
SGS  kinetic  energy  and  the  filtered  energy  dissipation  rate  for  the  SGS  velocity.  Therefore, 
LES  contains  the  essential  information  to  model  the  degree  of  nonequilibrium. 

The  results  also  show  that,  contrary  to  Kolmogorov’s  refined  hypotheses,  the  filtered  dissipa¬ 
tion  XL  and  ejr,  alone  are  not  sufficient  to  describe  the  SGS  scalar  and  velocity  because  the 
hypotheses  are  based  on  equilibrium  spectral  transfer. 

•  The  filtered  scalar  and  energy  dissipation  rates  cannot  be  modeled  accurately  using  the  spec¬ 
tral  equilibrium  assumption.  Nonequilibrium  effects  must  be  considered.  The  velocity  field 
affects  the  SGS  scalar  variance  and  filtered  scalar  dissipation  primarily  by  altering  the  degree 
of  spectral  nonequilibrium  and  the  SGS  time  scale  of  the  scalar,  respectively.  The  effects  can 
be  parameterized  by  the  SGS  kinetic  energy  and  the  filtered  energy  dissipation  rate.  There¬ 
fore,  the  models  for  the  SGS  scalar  variance  and  filtered  scalar  dissipation  can  be  further 
improved  by  including  the  effects  of  the  SGS  kinetic  energy  and  the  filtered  energy  dissipa¬ 
tion  rate. 
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In  collaboration  with  Professors  Pope,  Givi,  and  Pitsch  we  are  currently  working  on  testing 
models  and  LES  predictions  (FDF,  etc.)  and  are  exploring  the  possibilities  of  using  the  above 
suggestions  to  improve  models.  For  example,  the  measured  SGS  mixing  time  increases  with  the 
SGS  variance,  whereas  that  used  in  the  model  does  not.  This  model  behavior  does  not  allow  the 
SGS  variance  in  LES  to  fluctuate  to  the  large  values  observed  in  the  experiments.  In  addition,  the 
conditional  mixture  fraction  FMDF  obtained  in  LES  appears  to  have  the  bimodal  shape  at  SGS 
variance  values  much  larger  than  those  in  the  experiments.  These  model  behaviors  may  have  a 
strong  influence  on  the  predicted  reaction  zone  structure  and  will  be  addressed.  Further  activities 
in  this  area  will  be  included  in  future  project  reports. 
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Interactions  with  researchers  at  AFRL 

During  the  period  of  this  grant,  the  PI  initiated  contacts  with  several  researchers  at  AFRL  to 
seek  opportunities  for  collaborations.  Potential  topics  of  collaborations  were  identified  as  a  result  of 
the  subsequent  discussions  with  the  researchers.  In  particular,  the  discussions  with  Dr.  Campbell 
Carter  initiated  by  the  PI  in  2004  resulted  in  a  collaborative  project.  In  the  summer  of  2005,  the 
PI  spent  several  weeks  at  AFRL  as  part  of  the  collaboration. 

This  collaboration  focuses  primarily  on  issues  of  using  measurements  to  improve  of  large  eddy 
simulation  (LES)  of  nonpremixed  turbulent  combustion.  In  PDF-based  LES  approaches  the  subgrid- 
scale  mixing  of  multiple  scalars  must  be  modeled.  Current  mixing  models  are  based  primarily  on 
knowledge  gained  from  two-stream  mixing  problems  (e.g.,  fuel  mixing  with  oxidizer).  However, 
in  a  reacting  flow  at  least  three  scalars  are  involved  (the  third  is  a  product).  Therefore,  under¬ 
standing  of  three-stream  SGS  mixing  is  important  for  modeling  mixing  in  nonpremixed  turbulent 
combustion.  As  a  step  toward  understanding  SGS  mixing  of  multiple  reactive  scalars,  the  PI  and 
Dr.  Carter  will  study  the  SGS  mixing  in  a  three-stream  non-reacting  jet.  The  jet  nozzle  consists 
of  an  axisymmetric  jet  and  an  annulus  from  which  acetone-doped  air  and  ethylene  are  issued  into 
an  air  co-flow.  Laser  diagnostics  (planar  laser  induced  fluorescence  and  Rayleigh  scattering)  are 
employed  to  obtain  images  of  the  species.  This  study  will  provide  a  basis  for  future  investigations 
of  multi-scalar  SGS  mixing  in  turbulent  flames. 
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